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Factorizations  of  Normalized 
Totally  Positive  Systems 

J.  M.  Carnicer  and  E.  Mainar 


Abstract.  The  de  Casteljau  algorithm  for  evaluation  of  Bezier  curves 
can  be  generalized  to  curves  generated  by  any  normalized  totally  positive 
basis.  The  construction  of  this  algorithm  is  based  upon  a factorization 
of  the  system  as  a product  of  bidiagonal  stochastic  matrices  of  functions. 
These  factorizations  depend  on  a selection  of  a sequence  of  rectangular 
bidiagonal  matrices  of  decreasing  dimensions. 


§1.  Introduction 

The  Bernstein  basis  6"(f)  :=  (”)(1  — can  be  used  for  defining  a Bezier 

curve 

7(t  ):=£>TO,  *€[0,1]. 

i=0 

By  means  of  the  degree  raising  technique,  we  can  express  the  Bezier  curve  in 
terms  of  the  Bernstein  basis  of  one  higher  degree:  j(t)  — Qib™+1(t), 

t € [0, 1].  Indeed,  the  relations 

TO  = ^Tb?+1W  + ^TT6"+T  W-  < = 0, . . . ,n,  (1.1) 

can  be  written  in  matrix  form  as 


(TO-,*s)  = ( bn0+1,...xt\)A, 


(1.2) 


where  A is  an  (ra  + 2)  x (n  + 1)  nonnegative  stochastic  bidiagonal  matrix.  Such 
a matrix  can  be  written  as: 


(l  0 

| ai  1 — ai 


A = 
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a2 


° 
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0 an 
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0 < a;  < 1, 


(1.3) 
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Equality  (1.2)  corresponds  to  the  choice 

«»  := — rr>  i = (1.4) 

n H- 1 

Using  (1.2),  we  can  write 

7 (t)  = (*S , Pn)T  = (6S+1, . . • , 0^1(Po,  • ■ • , Pnf, 

which  proves  that  the  new  control  polygon  is  given  by 

(Qo,...,Qn+l)T:=^(Po,.-.!Pn)T.  (1.5) 

On  the  other  hand,  the  de  Casteljau  algorithm  for  the  pointwise  evalua- 
tion of  the  curve  is  based  on  the  following  well-known  recurrence  relations 

6"+1(t)  = Ai_1(t)6r_1(t)  + (l-A<(t))6"(t),  j = 0, . . . ,n  + 1,  (1.6) 

where  A i(t)  :=  t for  i = 0, . . . , n,  A_i (f)  :=  0,  and  An+i(t)  :=  1.  Indeed,  we 
can  write  (1.6)  as 

= (i.7) 

where  A(f)  denotes  the  nonnegative  stochastic  bidiagonal  matrix 
/l-Ao(t)  Ao(t)  \ 

A (t)=  ••.  • (1-8) 

\ 1-An(<)  A n{t)J 

Then,  starting  with  a Bezier  curve  7(f)  = Qtb™+]  (t),  (1.7)  gives 

7(f)  = (f>0+1.  • • ■ , OWo,  . . . , Qn+l)T  = (65,...,  6K)(Po(0,  • ■ • - 

where 

(Po(f),  • • • , Pn{t))T  :=  A(t)(Qo,  • • . , Qn+ if.  (1.9) 

Equality  (1.9)  describes  the  first  step  of  the  de  Casteljau  algorithm  for  the 
evaluation  of  7 (t). 

Bernstein  bases  are  totally  positive  on  [0, 1].  In  this  paper  we  shall  prove 
that  similar  properties  hold  for  any  totally  positive  basis  of  functions.  Let  us 
recall  that  a totally  positive  matrix  is  a matrix  such  that  all  of  its  minors  are 
nonnegative.  A totally  positive  system  of  functions  defined  on  / is  a system 
(uq,  . . . , un)  such  that  ('Uj(ti))o<i,y<n  is  totally  positive  for  all  <0  < • • ■ < tn  in 
I.  A normalized  totally  positive  (NTP)  basis  (uo, . . . ,un)  is  a totally  positive 
system  of  linearly  independent  functions  such  that  )T)?=o'u»  = L 

Given  an  NTP  basis  (ug+1 , . . . , u^+i)  on  an  interval  I,  and  a nonnegative 
stochastic  (n  + 2)  x (n  + 1)  matrix  A of  rank  n + 1,  we  shall  consider  the  system 
of  functions  defined  by 


«,...,<)  :=K+1,...,<+i)A 


(1.10) 
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Starting  from  a curve  7 (t)  = Plurf(t))  t e I,  clearly  (1.10)  allows  us 
to  express  it  as  7 (t)  = YaIq  <3i«"+1(t),  where  the  points  Qj  are  defined  by 
(1.5).  In  this  paper  we  will  derive  from  (1.10)  the  existence  of  nondecreasing 
functions  Ao, . . • , A„  with  values  in  [0, 1]  such  that 

K+1(f),...Xil(*))  = W(t),...,<W)A(f),  tel.  (1.11) 

The  matrix  Aft)  is  defined  from  A,  as  in  (1.8).  Starting  with  a curve  7 (t)  — 
Sfco1  Qi<+\t)>  t e I,  v?e  will  be  able  to  write  it  as  7 (t)  = J27=o  Pi(t)u?(t), 
where  the  points  Po(t), . . . , Pn(t)  are  again  given  by  (1.9).  On  the  other  hand, 
we  shall  check  that  (1.10)  implies  that  (uq,  . . . , u")  is  an  NTP  basis  on  I.  It 
will  therefore  be  possible  to  iterate  this  process.  Doing  so,  we  shall  obtain  a 
de  Casteljau  type  algorithm  for  the  evaluation  of  7 (t). 

Pottmann  and  Mazure  in  [5,6,7]  developed  generalizations  (1.11)  of  the 
de  Casteljau  algorithm  for  Tchebycheffian  curves.  Here  we  show  that  these 
generalizations  can  be  also  obtained  for  any  curve  generated  by  an  NTP  basis. 

We  observe  that  for  each  value  of  t,  the  point  Pf  (t)  is  a convex  combina- 
tion of  two  consecutive  points,  obtained  in  the  previous  step  of  the  algorithm. 
Therefore,  these  algorithms  can  be  seen  as  corner  cutting  algorithms  for  curve 
evaluation  [4]. 


§2.  Recurrence  Relations  for  NTP  Systems 

The  following  proposition  allows  us  to  describe  the  generalization  of  formulae 
(1.6)  to  any  NTP  systems  related  by  a matrix  (1.3).  First  we  need  to  show 
the  following  auxiliary  result. 

Lemma  2.1.  Let  («q+1,  . . • , un+i)  he  an  NTP  basis  of  functions  defined  on 
I and 

«>•••»<)  :=(«o+1.*".«S+iM>  (2-1) 

where  A € ]R,(n+2)x(n+1)  jg  of  the  form  (1.3).  Let  Ci  :=  {t  e I | u"(t)  ^ 0}. 
Then 

(i)  u"+i{t)/ui(t),  t € Ci,  is  a nondecreasing  function, 

(ii)  ai+1u^(t)/uf(t)  e [0, 1],  for  all  t e (7*. 

Proof: 

(i)  Since  A is  bidiagonal,  we  can  write 


<(t)  = (l-o0<+1(t)  + ai+i<+1(t),  tel,  i = 0, . . . ,n.  (2.2) 

Observe  that,  since  A is  nonnegative  and  (uj+1, . . . , «"+])  is  NTP,  then 
u"(f)  > 0,  for  all  t e Ci . Moreover, 


«"(«) 


<++ixW 


= (1  - Oj) 


<+1(f) 


Ci(s) 


>0, 


t < s, 


(2.3) 


because  (wq+1,  . . . , «"£])  is  totally  positive.  Formula  (2.3)  implies  that 
/uf  is  nondecreasing  in  Ci. 
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(ii)  Using  (2.2),  we  can  write 
«i+i  (<)  _ 


0 < 


«i+i«?+i  (*) 


«?(*) 


, - < 1,  teCi.  □ 

(l-aiK+1(0  + ai+1i*?+1(t)  - 

The  following  proposition  is  devoted  to  showing  that  formula  (1.11)  holds 
for  NTP  bases. 


Proposition  2.2.  Let  («q+1,  . . . , u”+|)  and  (ng , . . . , u")  be  two  NTP  bases 
of  functions  on  I related  by  (2.1),  where  A is  a matrix  (1.3)  (rank  A = n + l). 
Let  Ci  :=  {t  £ I \ n"(t)  ^ 0}.  Then  the  functions  A*  : I — ► 1R,  i = 0, . . . , n, 
defined  by 

x.rt\  f «i+i inf  {“"+11  («)/«“(«)  1 s € Cj,  if  <(s)  = 0,  Vs  < t, 

1 ' \ Oj+i  sup  {«"+i1(s)/n"(s)  j s £ Ci,  s < f},  otherwise, 

(2-4) 

are  nondecreasing,  and  satisfy 

0 < Aj(t)  < 1,  Vf  € I,  t = 0, . . . ,n.  (2.5) 

Furthermore,  if  we  use  definition  (1.8),  then  (1.11)  holds. 


Proof:  Since  (no,...,u")  are  linearly  independent,  then  C\  ^ 0 for  all  i. 
Therefore,  by  Lemma  2.1  (ii),  we  can  define 

«ec<|e[0,i]. 

If  the  condition  «"(s)  = 0,  Vs  < t,  does  not  hold,  then  the  set  {s  6 C,  | s < t } 
is  nonempty  and  by  Lemma  2.1  (ii),  we  can  define 

f un+1(s ) 1 

Ki  < A i(t)  :=  ai+i  sup  | s £ Ci,s  < < 1. 

We  have  seen  that  A*,  * = 0, . . . ,n,  are  well-defined  and  that  (2.5)  holds.  In 
order  to  see  that  A*(t)  are  nondecreasing,  let  us  observe  first  that  if  {s  € Ci  | 
s < tj}  = 0 and  {s  £ Ci  \ s < <2}  7^  0,  then  <i  must  be  less  than  <2-  Therefore, 
we  only  have  to  show  that  A;(ti)  < A;(<2)  only  for  all  ti  < t%  such  that  there 
exists  some  s <t\  with  u*(s ) ^ 0.  We  observe  that  {s  £ Ci  | s < t\}  C {s  £ 
Ci  | s < <2}.  Therefore,  by  Lemma  2.1  (i), 


ai+1  sup 


«?(«) 

We  now  establish  the  relation 


ai+i«r+i(<)  = Aj(t)ii"(t),  Vt  € I,  i = 0,...,n.  (2.6) 

If  u"(f)  = 0,  then  by  (2.2),  ai+iu^(  (t)  = 0,  and  (2.6)  trivially  holds.  Other- 
wise, we  have  t £ {s  £ Ci  \ s < t},  and  by  Lemma  2.1  (i), 

Ai(f)  = «i+i<+i1(f)/«?(t), 

so  (2.6)  is  confirmed  again. 

Finally,  using  (2.6)  and  (2.2)  we  can  write 

Ai-r (*)«?_, (t)  + (1  - A i(f)K(f)  = aiu?+l(t)  + «?(«)  - a>+1  <+/(<)  = <+1(f) 
for  all  t £ I and  i — 0, . . . , n.  □ 
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§3.  A Generalization  of  the  de  Casteljau  Algorithm  for  NTP  Bases 

Given  an  NTP  basis  (uq,.  ..,«”)  of  a space  Un  of  functions  defined  on  I, 
we  can  obtain  a sequence  of  NTP  bases  ( Uq , . . . ,u%)  of  (k  + l)-dimensional 
subspaces  Uk  by  the  recurrence 

(ukQ(t),...,ukk(t))  :=  (uJ+1(i)1...,«J+}(t))Alfc+i,  k = n—l,n—2, . . . ,0,  (3.1) 

where  Ak+i  € ]R(fc+2)x(fe+1)  is  a matrix  of  type  (1.3),  rank  Ak+i  = k + 1. 

In  fact,  since  Ak+i  are  nonnegative  bidiagonal  matrices,  it  easily  fol- 
lows, using  Theorem  2.3  of  [1],  that  Ak+\  is  totally  positive  and,  using  the 
Cauchy-Binet  formula  (formula  (1.23)  of  [1]),  that  the  systems  (3.1)  are  totally 
positive.  Taking  into  account  that  (uq+1,  . . . , u^_\)  is  normalized  and  Ak+\  is 
stochastic,  we  derive  that  the  systems  (3.1)  are  also  normalized.  Furthermore, 
formula  (3.1)  relates  two  bases  if  and  only  if  rank  Ak+i  = k + 1.  Observe  that 
rank  Afc+i  < k + 1 if  and  only  if  there  exist  1 < i < j < k such  that  a*+1  = 1 
and  ak+1  — 0. 

Let  us  observe  that  the  subspaces  Uk  form  a chain,  that  is, 

Un  D Un~x  3 ■ ■ O W1  D M°  = span{l}. 

Moreover,  since  («[])  is  an  NTP  basis  of  14°  then  f°r  ah  t € I- 

By  Proposition  2.2,  the  bases  of  (3.1)  are  related  by 

Kfc+1(t),...,^+l(t))  = (4(t),...,4(t))Afe+1(t),  t€l,  (3.2) 

where  Ak+i(t)  is  a matrix  of  type  (1.8).  We  shall  denote  by  A*+1(f)  the 
(i  + l,i  + 2)  entry  of  Ak+i(t).  The  recurrences  (3.1)  and  (3.2)  give 

(«o(t), . . .,uk(t))  = («o(t),  • • • ,<(f))An  • • • Ak+2Ak+i,  tel,  (3.3) 

and 

(«o(0,-*-,«fc(0)  = Ai(t)A2  (*)•••  Afc(t),  tel,  (3.4) 

for  k = 0, . . . , n,  with  the  convention  An  ■ ■ ■ Ak+ 1 equals  the  identity  matrix 
when  k = n and  Ai(t)  • • • Ak(t)  equals  the  scalar  constant  1 when  k = 0. 

Formulae  (3.4)  can  be  interpreted  as  a factorization  of  the  NTP  system 
( uk , ...,«£)  as  a product  of  bidiagonal  stochastic  matrices  of  functions. 

Let  us  summarize  all  the  conclusions  in  the  following  theorem. 

Theorem  3.1.  Let  (iifi, . . . ,«")  be  an  NTP  basis  of  functions  defined  on  I. 
Let  Ak  e JR{k+1)xk , k = 1, . . . ,n,  be  matrices  (3.1)  of  maximal  rank.  Define 
NTP  systems  (uk, . . . , u£),  k = 0, . . . ,n  — 1,  by  (3.1)  (or  equivalently  by  (3.3)). 
Then  there  exist  matrices  A k(t)  of  type  (1.8)  whose  (i  + l,i  + 2)  entry  X k(t) 
is  nondecreasing  on  I and  with  values  in  [0,1],  k = 1, . . . , n,  such  that  (3.2) 
and  (3.4)  hold.  In  particular, 


K(t),...,<(t))  = A1(t)...An(f),  'it  el. 


(3.5) 


6 


J.  M.  Camicer  and  E.  Mainar 


Moreover,  for  any  control  polygon  Pq  ■ ■ ■ Pn  consider  the  following  generaliza- 
tion of  the  de  Casteljau  algorithm: 


for  j = 0,1,..., n 


m)  ■■=  pj 

for  i = n — 1, . . . , 1,0 
for  j = 0,1,..., i 


Pj(t)  :=  (1  - A j+1(i))P;+1(l)  + A*+1(t)P*+}(t) 


At  each  step  we  have 


l(t)  = ^PJ(t)u}(t),  tel,  i = 0, . . . ,n.  (3.6) 

j=o 


In  particular  7 (t)  = Pq(1)  for  all  tel,  that  is,  this  generalized  de  Casteljau 
algorithm  reconstructs  the  curve  from  its  control  polygon. 

Proof:  The  existence  of  the  matrices  A k(t)  of  type  (1.8),  satisfying  (3.1)  and 
(3.2)  follows  from  Proposition  2.2.  From  the  algorithm  we  see  that 

(m  • • pmT  = a i+1(t)(p’+i(<), ....  pi$(t))T, 


and  by  (3.5)  we  can  write 


7(t)=(uS(t),...,K(t))(p0,...,pn)T  = 

Ai(t)  ■ ■ ■ Ak(t)Ak+i(t)  ■ ■ ■ A n(t)(P0, ....  Pn)T  = 
(ui(t),...,ukk(t))(p^(t),...,pkk(t))T.  □ 


Example  3.2.  When  applying  Proposition  2.2  to  (1.1)  or  (1.2),  the  func- 
tions that  we  obtain  are  A ,(t)  = t,  i = 0, . . . ,n.  Hence  we  obtain  (1.7),  and 
the  corresponding  algorithm  described  in  Theorem  3.1  is  just  the  classical  de 
Casteljau  algorithm  for  polynomials.  Of  course,  any  other  choice  of  a sequence 
(ylfc),  k = 1, . . . , n,  of  nonnegative  stochastic  matrices  of  maximal  rank  could 
lead  to  another  de  Casteljau  type  algorithm.  For  instance,  if  we  consider  the 
Bernstein  basis  (6q,  , 62)  °f  degree  2,  the  matrix 


defines  aNTP  basis  ((1  — 1)(3— 1)/3,  t(4  — 1)/3)  on  [0, 1].  This  system  generates 
a subspace  of  quadratic  functions,  different  from  the  polynomials  of  degree  less 
than  or  equal  to  1.  Furthermore,  the  functional  matrices  obtained  by  applying 
Proposition  2.2 


/ 3(1— t) 
[ 3-f 

V 0 


2 1 

3 -t 


4(1-0 

4— t 


0 

3 1 
4— t 


Ai(t)  = ( ^),  A2(t) 
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lead  to  a corner  cutting  algorithm  different  from  the  classical  de  Casteljau 
algorithm. 

In  [2],  it  was  shown  that,  in  any  space  with  an  NTP  basis,  there  exists 
a particular  NTP  basis  called  the  normalized  B-basis  which  has  the  optimal 
shape  preserving  properties  among  all  NTP  bases  of  the  space.  In  Theorem 
4.3  of  [3]  it  was  shown  that  if  (mq+1,  • • • >un+i)  's  a normalized  B-basis  of 
an  (n  + 2)-dimensional  space  and  (ng,.  . . , u")  is  a B-basis  of  an  ( n + 1)- 
dimensional  subspace,  then  there  exists  a matrix  A (1.3)  such  that 


Thus,  B-bases  provide  good  examples  of  when  Theorem  3.1  can  be  applied. 
In  the  case  of  polynomial  spline  spaces  (see  [2]),  the  normalized  B-basis  is 
precisely  the  B-spline  basis. 

Let  T = {t0  = •••  = tk—l  < tk  — ' — tn  < tn+i  = * * * — tn_|_fc{, 

tj  < tj+fc,  for  all  i,  be  an  extended  knot  sequence  and 

~ i^i+k  ~ ti)[ti>  • • • , ^i+fc]("  — ^)-f  , t € [^Oi  ^n+l]i  * — 0,  . . . , H, 

the  associated  B-spline  basis  of  the  space  <S£.  Let  us  insert  a knot  r in  T such 
that  tj  < t < tj+i  (if  r = tj  then  the  multiplicity  of  tj  must  be  less  than  k) 
and  define  a new  sequence  of  knots  T 

(U,  0 <i<j, 

tj  :=  < r,  i = j + 1, 

t tj_i,  j + 2<*<n  + A:-|-l. 


The  normalized  B-bases  of  Slf  and  Sff  are  related  by  a matrix  (1.3)  with 


r o,  o < i < j — k + 1, 

on  :=  (tj+fc-i  - r)/(ti+k-i  - tj),  j - k + 2 < i < j, 

L 1,  j + 1 < i < n. 

Applying  Proposition  2.2  to  both  B-spline  bases,  a relation  (1.11)  is  obtained. 

In  order  to  obtain  a generalized  de  Casteljau  algorithm,  we  first  remove  suc- 
cesively  all  interior  knots  until  we  arrive  at  the  Bernstein  basis.  Then  we  can 
continue  with  the  steps  of  an  evaluation  algorithm  for  polynomials  (e.g.,  the 
de  Casteljau  algorithm).  We  illustrate  this  procedure  with  a simple  example: 

Example  3.3.  Take  T :=  (0, 0, 1/2, 1, 1),  T :=  (0,0, 1,1).  The  associated 
B-spline  bases  are  related  by  a matrix  (1.3): 


t G [0,1]. 


Using  Proposition  2.2,  we  obtain  that  these  bases  are  also  related  by  (1.11) 
with  a matrix  (1.8),  where 

A0(t)  = min(l,t/(l  - f)),  Ai(t)  = max(0,  (2f  - 1) /t),  t € [0, 1]. 
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The  evaluation  algorithm  for  7 (f)  :=  Ya=o  P,N'^{t)  can  be  described  as 
follows.  First  compute 

Po(t)  ~ (1  - Mt))P0  + Xo(t)Pu  Pl(t)  :=  (1  - A 1(t))P1  + Xi(t)P2, 

and  then  7 (t)  = (1  — <)P01(t)  + tPy(t).  Note  that  the  last  step  of  the  algo- 
rithm corresponds  to  the  de  Casteljau  algorithm.  Of  course,  this  algorithm  is 
different  from  the  classical  de  Boor-Cox  algorithm  for  evaluation  of  B-splines. 

Acknowledgments.  This  research  has  been  partially  supported  by  Research 
Grant  DGES  PB96-0730.  We  thank  the  referee  for  his  valuable  suggestions. 
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